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Abstract 

A diagonal Bianchi Type III space-time is treated, both at the classical and quantum 
level, in the context of Hofava - Lifshitz gravity. The system of the classical equations of 
motion is reduced to one independent Abel's equation of the first kind. Closed form solution 
are presented for various values of the coupling constants appearing in the action. Due to 
the method used, solutions of Euclidean, Lorentzian and neutral signature are attained. The 
solutions corresponding to A 7^ 1 are seen to develop curvature singularities as the other 
constants approach their Einsteinian values, in contrast to those with A = 1 which tend to 
the known Einstein gravity solutions. At the quantum level, the resulting Wheeler-DeWitt 
equation is explicitly solved for A = 1, cr = and A = ^. The ensuing wave- functions diverge 
in the Einsteinian limit. 
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1 Introduction 



It could be stated that, an important feature of general relativity is the covariance of the theory 
under four dimensional diffeomorphisms of the spacetime manifold. Recently, a new theory of 
gravity was proposed, with the aim of being complete in the UV [I]. Its basic assumption is 
the admittance that the aforementioned covariance is not a fundamental property of the theory, 
but arises rather accidentally in the context of lower energies. The resulting geometry is then 
exhibiting a non isotropic time and space scaling invariance 

As a consequence, one must start from an action functional and equations of motion which involve 
higher derivatives in the spatial coordinates. On the other hand, it is an advantage that there are 
no higher derivatives in time. A theory was constructed for various values of the critical exponent 
z (2 = 2 [2], 2 = 3 and 2 = 4 

In the context of this theory one has to consider a four dimensional different iable manifold, M., 
with a codimension one foliation T and admit as the "gauge invariance" group the foliation 
preserving diffeomorphisms, i.e. coordinate transformations of the restricted form 

t = f{t) , 5^=g\t,x^). (1.1) 

With this reasoning the author of [TJ modified the 3+1 decomposed Einstein - Hilbert action by 
adding one extra coupling constant in the kinetic term that destroys four dimensional covariance; 
As far as the potential part is concerned, extra terms were added under the condition of detailed 
balance. The latter results in a potential term of the lagrangian density, which up to a coupling 
constant is 

r rp f~iijkl rp 

^potential ^ij^ J-^kl 

with G'*-?'^' the generalized Wheeler - DeWitt supermetric and Eij a three tensor that is obtained 
by a variational principle ^/gE^^ = ^^J-^'''-^ for some action W pQ. 

til J 

There are many papers studying certain cosmological implications of the theory, mostly for a FRW 
metric, either free of matter or with a scalar field (see for example or for an excellent review 

on Hofava - Lifshitz cosmology, [10] and the references therein). In this paper we are interested 
in following the point of view of this theory to obtain solutions in the classical context, for a 
diagonal Bianchi type III model in vacuum, which, as we know, has an anisotropic spatial metric, 
and compare with analogous solutions from Einstein's equations. Moreover, we proceed with the 
canonical quantization of the axisymmetric case and derive the Wheeler - DeWitt equation for 
this cosmological model, giving a solution for specific values of some of the coupling constants. 
The results obtained are also compared to those obtained by quantizing the Einsteinian action. 
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2 Equations of motion 



Our starting point is the action [3] 

S = J dt£x^N (^a{KijK'^ - XK^) + -C^.^C^'^ + -f^RuR\.p + CR)R{ + iiR^ + iR + (t^(2.1) 

where g is the determinant of the three dimensional metric Qij, N is the lapse function, Kij = 
jj^{gij — Ni-j — Nj-i) the extrinsic curvature corresponding to the spatial metric, K its trace, 
(jpq _ ^pk<-(^]^'i^ — ^R5\)-k the Cotton - York tensor density and Rij, R the Ricci tensor and scalar 
respectively. Each term in f l2.ip comes with its own coupling constant a. A, /3, 7, ^5 i and a. 
The two parts of the kinetic term are distinguished by the existence of A. For A = 1 we get the 
kinetic term of Einstein's theory and the difference is restricted to the extra terms in the potential 
part and the different coupling constants that bind them together. We note here that ; stands for 
covariant differentiation with respect to the spatial metric and all the latin indices run from 1 to 
3. 

Variation of the action by 5gmn results in the following spatial equations of motion 

o nr qkl qkl 

■^9 '^9 ^9 



\ksl 

^qkm/^ n ri i ^nkn/-~i m t) \ \ 
e K-k + e '^q ti;k) + 



'^9\ 

Jpk ( 

-^^NiRrRl, + R^RIp) + ^(iVi?,fc;p); J e"^^'^ + e'r>kgmqgsn _ 

_^lpkgsqgran _ ^mpk ^nq ^sl ^ ^ _ 



7 / N^iplgsnj^rn ^ ^ ^^pm gSn p,l _ ^ ^^ps gVan _ ^ ^^pn gSl ^rn ^ ^ ^) 



\psl 



+ Ne''^''RlR\ - Ne'^'^Rl^R^ + (m o n) 1 + 

+ ^R'^Rijg"''' - 2CNR']R^'' + 

+^R^g'^'' - 2r]NRR"''' + 7]{N R) ,^^0^'^'^'' + 

+ ^g^n ^ ^Rgmn _ ^NJ^rnn ^ ^pq^^^Qpqmn ^ q ^2.2) 

Where (m n) stands for the repetition of the terms in the parenthesis followed by an interchange 
in the m, n indices and G^'^"^"' = gP"^gi"- -\- gP^-gi"^ — 2gP'^g"^"' is the covariant supermetric. Variation 
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by 5Ni yields the linear constraints 



(2.3) 



and finally variation by the lapse function N leads to the quadratic constraint 



E"^ = -aiKijK'^ - XK') + /SdjC' + -f—RuRkj + CR^JR'' + riR' + iR+a = Q (2.4) 

V 9 



3 Diagonal Bianchi type III, classical case 

As it is known [TT] group automorphisms can be used for simplifying Einstein's equations, since 
the automorphism generators are Lie point symmetries of the system. In our case we will consider 
a diagonal scale factor matrix, which, by the procedure described in [T2] and [13] , must be written 
as 

/ gni(t)+«3(t) Q g \ 

g g«2(t)+-it3(t) g 

e"3(*) 



7a/3 



V 
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In this particular parametrization of the 7a^'s the automorphisms generators are cast into canon- 
ical form, so that the ensuing equations become of first order with respect to the corresponding 
velocities. The invariant basis 1-forms are 





^0 e-^ 









1 








given by 





/g"3{t) 



9t3 



a P 



\ 










-2x+ui{t)+U3{t) 








^U2{t)+U3{t) 



One can easily check that this metric corresponds to a three dimensional conformally flat space. 
This means that the corresponding Cotton - York tensor is zero. Moreover, one can also see that 
the covariant derivative of the Rij tensor is also zero. Since the invoked coordinate transformations 



in [TT] are precisely of the form fll.ip we can, without any loss of generality, take the shift vector 
A^* to be zero. Then the linear constraint (12. 3p yields the equation 



ae 



N{t) 







(3.1) 



while the quadratic constraint (12. 4p becomes 

4N{ty (2C + 4:7] + e"«(*)(-2e + ae"^^*))) + ae'^^^ ((A - l)(Mi(t)' + ti2(t)') + 
+ (3A - I){2u2it)u3{t) + 3u3itf + 2ui{t)u3{t)) + 2Xui{t)u2{t)) = 0. 



(3.2) 
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Finally, the three nonzero spatial equations reduce to 

E]]j = -4N^{2C + At]- ae^"^) + ae^"^ (^Nu^ + N{ul + ul - 3ul - Au^)^ + 
+aAe2"3 (n{ui + U2 + 3%)^ + 4(-iV(Mi + U2 + 8^3) + N{u^ + U2 + = (3.3) 



Efjj = -4N%2C + 4r]- (xe^"-') - iae^'^'N ((A - l)ui + \u2 + (3A - 1)%) + 
+aNe^'''' ((A - l)ul + (A + l)ul + 6XU2U3 + 3(3A - 1)^^ + 2(A - 1)mi(m2 + 3^3) + 

+4(A - l)ui + 4Am2 + 4(3A - l)ii3) = (3.4) 

Effj = (2C + 477 + e"'(-2e + fxe"^)) - 4«iVe2"^ (Ami + (A - 1)^2 + (3A - 1)^3) + 
aA^e^"^ ((A + l)uj + (A - l)ul + 6(A - 1)m2M3 + 3(3A - l)ul + 2mi((A - 1)^2 + 3Am3) + 

4Aiii+4(A-l)M2 + 4(3A-l)M3) = (3.5) 

It is useful to observe that in equations 03.31) . fl3.4p and f l3.5p there is no actual need to take the 
square root of N{t)'^: dividing these equations by N(t) we see that only N(tY and = 2^^(^^) 
appear. Therefore if we do not "gauge" fix the lapse a priori, solutions with any signature will be 
attained. Equation (13. ip implies that ui{t) must be a constant. Upon setting ui(t) = ci, where ci 
is a real constant, equations Ej]j = and E]]j = become identical. At this point we can solve 
equation (13. 2p with respect to the lapse function and substitute its value into the independent 
spatial equations (13. 3 p and (13. 5p . This action leads to one final independent equation in terms of 
the remaining functions U2{t) and u^lt). Since this equation is of first order with respect to U2{t) 
(i.e. U2{t) does not enter it), while of second order in u^lt), it is almost mandatory to choose u^lt) 
as time. With the specific choice u^lt) = ln(t) and upon setting U2{t) = J uj{t)dt we arrive at a 
final equation of first order in u: 

._ _ t(A-l)(2AC + 4Ar/ + t((Tt-(A + l)0) 3 _ 4r/(3A - 1) + 3C(3A - 1) + t(2at - ^ - 3A0 2 
^~ ~ 2(3A- l)(2C + 4r7 + t(at-20) ^ 2{2C + At] + t{at - 2^)) ^ 

18AC + 36A?7 + t{5at - 5g - QAQ ^^ 3(3A - 1)(2C + 4r] - ^t) 

This is an Abel equation of the first kind. We did not manage to identify the above equation with 
any known integrable class, so in order to simplify the situation, in the following subsection we 
will set A = 1. 

But, before we do that, we have to investigate the situation where the coefficient of N{ty in (13. 2p 
vanishes, so that this equation does not determine N{t). This happens for 



usit) = ln( ^ ) (3.7) 

(J 

and then the quadratic constraint (13. 2p dictates 

(A-l)M2(t) = (3.8) 
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Since the case A = 1 will be treated in the next subsection, we will here restrict ourselves to the 
case U2{t) = C2 with C2 a real constant. With U2 being constant and by the help of (13.71) . equation 
()3.5p is identically satisfied, while the other independent spatial equation (13. 3|) (= (13. 4p ) becomes 



2(( + 2,)-MVlEM±M)!j,VW»=0 (3.9) 

Since the lapse function cannot be zero, in order for a solution to exist, we have to reduce the 
number of independent coupling constants. We solve (13. 9p for a and find that a = 2{(^^+2r]) both 
values of u^. The spatial metric has now become 



/2(^ \ 



9ij 



2^^e^i-2- 



2^^e^2 

and one can write the line element as 



V 



ds^ = ede + 2^^^j^dx' + e-^^dy'' + dz" (3.10) 

where e = ±1 and the following simplifications have been made: since the lapse has not been 
defined by (13.21) . a reparametrization of time has been used to set N{t) = e and the spatial 
coordinates {y, z) have been rescaled in order to absorb the constant factor K:^^) together with 
the non essential constants ci and C2- As we see, only the coupling constants ^, Q and r] appear 
in (I3.10p . This is expected by the assumption (13. 7p and its consequences; since the ensuing 
pseudometric is static A is excluded from it. Moreover, this is a somewhat isolated solution, in 
the sense that in the Einsteinian limit (C — )■ 0, 77 — )■ 0, ^ — )■ a) it develops a curvature singularity, 
because the Ricci scalar is i? = This particular space admits six killing vectors 

^ 9 d d d 

oy OZ ox oy 

The first three are what we would expect from the isometry group and satisfy C^^a"" = 0, with 
A = 1,2,3. The fourth depends on the metric components and is of the form investigated in [15] 
and [16]. Usually this kind of vector does not exist as spacetime killing field, but is limited on 
the spacelike hyper-surface [13]. The non zero structure constants of the algebra closed by these 



six kilhng fields are Cl, = -C^, = Cf, = -Cl = Cl, = -Cl, = Cl, = -Cj, = Cj, = -Cj, = 
1. Finally, the metric exhibits the property that its Riemmann tensor has vanishing covariant 
derivative Rijkl\m = [E] (I, J., run from to 3 while | stands for 4d covariant differentiation). 
Therefore there are not higher derivative curvatures. 

3.1 Case A = 1 

Under the value A = 1, as already mentioned, the kinetic term of the theory is identical to the 
corresponding of Einstein's gravity. By this assumption, equation (13. 6p . conveniently becomes a 
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generalized Riccati equation. 

18C + 3677 - 14^t + bat 



-u 



3(2C + 477-et) 



2t(2C + 4:1] + t{af^ - 2^) t'^{2C + Ar] + t{at - 2^) 



whose solution is 



u{t) 



3(6c2et + 12c2(C + 2?7) + x/t) 



t{3c2at^ - 18c2(^t + C + 2r])- 2Vi) 
with C2 being the constant of integration. The lapse function becomes 

N{t) = 



{3.11] 



(3.12) 



ac2 



-Qc2at^ + 36c2(et + C + 2r]) + 4^/i 



(3.13) 



while the three dimensional metric is 



9tj 



ft 





\ 










^0 


3c2at2 - l8c2(et + C + 2r/) 


- 2yft 


t 


1 



From the above metric, we see that Ci is not an essential constant, since it can be extinguished by 
a simple scaling of the y coordinate. The corresponding line element is 



2P{t) 



t 



(3.14) 



where P{t) = 3c2at^ — 18c2{C,t + C + 2?]) — 2\/i. The integration constant C2 and the positive 
or negative value of P(t) in some open interval of t G M"*" determines the signature of the line 
element: (a) C2 > and P(t) > the signature is Euchdean, (b) C2 < and P(t) > the signature 
is Lorentzian and (c) C2 > and P{t) < the signature is neutral (— , +, +, — ). It is interesting 
that despite the non linearity of the equations, the above line element can be identified with the 
known solutions of Einstein gravity: by setting ( = rj = a = and ^ = a we arrive to vacuum 
diagonal type III ([E], [I3]), while setting ^ = 77 = and ^ = a we obtain the cosmological 
constant type III solutions ([IS])- 

As we mentioned before, we also have to study the case in which fl3.7p holds true. In view of 
(13. 7p . (13. 8p the quadratic constraint and (13. 5p become identically zero and we are left with the 
independent spatial equation 



a 



- ^{C + 2??) + we - MC + 2r?)) [n{2u2 + ill) - 2NU2) + 
+2a - MC + 2??) + We-MC + '^v)) = 



(3.15) 



We have the freedom to choose U2{t) = t and thus the above equation can be integrated for the 
dependent variable N{t). The result is 



N{ty 



1 - /i2e*+^i^2 



(3.16) 
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where 



= |i±^?EM±M'| (3.17) 



/U2 



a 

4^ 



i(2(C + 2r/)-^ia) (3.18) 



2a 

and C2 is the constant of integration. With these substitutions the spatial metric assumes the form 



9ij 



\ 

T^flTe^i-^^' 
V 77ZTe*y 



and the hne element can be written as 



ds^ = -dt^ + JJUdx' + e-'^'^dy'' + e'dz' (3.19) 

1 - /i2e* 

Again, depending on the sign of 1 — ^2^^ the solution can be both of Lorentzian or Euclidean 
signature: if ^2 is negative, then the metric is Lorentzian for all t G M, whereas if /i2 > the 
solution is Lorentzian for t < — In ^2 and Euclidean for t > — ln/i2. In the Einstenian limit — )■ a, 
C — )■ and ?7 — 7- (we don't set a — because it obviously leads to singularity), we are led to a 
singular pseudometic if we choose the minus solution in fl3.17p . On the contrary, for the plus value 
in f l3.17p we are left with a pseudometric that is solution of Einstein's equations plus a cosmological 
constant with A = The space-time with line element f l3.19p admits six killing vectors 

d d d d 

^1 = « = S3 = ^ — \- y-^ 

ay oz ox oy 
d z ^- d 



^6 = e2 Ve-* - //2^ - 62 v^e^* - //2 



The non zero structure constants of the corresponding algebra are: = — C|i = Cf^^ = — C|i = 
CI4 = -Cts = 1, CI, = -CI, = CI, = -CI, = i, Cfe = -CI, = ti2- Again the covariant 
derivative of the Riemmann tensor vanishes {Rijkl\m = 0). The constants /ii, /i2 appearing in the 



line element can both be seen to be essential by considering the curvature scalars R = — ^~^y^^ 
. Hgd. s„.e tKe sys.e. can L solved fo. .., tenns of these two cu^.e 
scalars, there is no spacetime coordinate transformation that can alter their values; thus they are 
essential also in the context of the present theory which allows less freedom i.e. transformations 

3.2 Case X = l 

Equation (13. 6p implies that we also have to distinguish the case A = |, since (3A — 1) appears 
in the denominator. This is the value for which Trace{Kj — XK6j) becomes zero. So, for that 
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specific value of A, we proceed in tlie same manner, and after solving (13.21) with respect to the 
lapse function and substituting in the spatial equations, we are left with the two independent 
equations produced by f l3.3|) and fl3.5|) 

(2(C + 2ri) + e"3(-4e + 2ae"^)) ul{u2 + M3) = (3.20) 
(2(C + 2r]) + e"«(-4^ + 2cre"«)) iilii^ = (3.21) 

When we subtract (I3.2ip from (I3.20p we end up with two possibilities. Either U2{t) is a constant or 
Us{t) is a constant whose value extinguishes the parentheses. It is easy to see that U2{t) = constant 
leads, upon setting X = ^, ui(t) = constant in (13.21) . to a zero lapse function. Thus we choose 



.^ , / 2f ± v/4f 2 - QaiC + 2n) \ 
Mt) = In I ^ j (3-22) 

Upon substitution all the spatial equations are satisfied and for the lapse function we have 

afe± v/4e-6(C + 2?7))' 

mr = — 7 — ^ ^ . Mt? (3.23) 

12a [-e + 6a(C + 2r/) t ^^4^ - 6a(C + 277) j 
with U2{t) an arbitrary function of time. 

By setting ui = '^^^^^^ 6o-(c+2??) _ , , . The spatial metric is 

written as 

f n n \ 



9ij 




















while, upon choosing U2{t) = t, the four dimensional line element is 

ds"^ = -vliy2dt^ + vidx^ + e'^'^dy^ + e^dz"^ (3.24) 

where again the necessary simplifications have been made. Once more, the relation of the values 
of the various coupling constants determines the signature via the sign of V2- The corresponding 
Riemmann tensor has vanishing covariant derivative and the two constants z^i, V2 are essential, 
since R = and RuklR^^^^ = ^^suff ■ ^^'^^^ ^^^^^ ^ ^ a, C ^ 0, 77 we are 

led to a singularity for the minus solution. On the contrary, the plus solution gives a Riemmanian 
metric which solves Einstein's equations plus a cosmological constant A = The killing vector 
admitted by (I3.24p are 

9 9 (9 d 

ay oz ox ay 

forming an algebra with the non zero structure constants C}^ = —C^i = Cf^^ = —Cfi = C|g = 

~^62 = ^34 = ~^t3 = ^56 = ~^65 = ^ ^^d Cfg = —C^2 = 2. 
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3.3 Case A^^l, ( = rj = a = 

Finally, we mention two choices for the parameters that lead to another solution. We leave A as 
an arbitrary constant different from unity and set C = r] = a = 0, or ( = —2// , a = 0. Both 
choices lead to the same differential equation for U2{t). If we adopt the choice U2{t) = J ^^^dt, we 
are left with 

\ 2 _ 1 

itu + — + (3A + 1)^2 + (1 + 9\)u + 3(3A - 1) = (3.25) 

3A — 1 

which has the solution 



u = F-\f{t)) (3.26) 



,,,1,-1 ( 3A-l + (A-l)/(t) \ _ 
4(A-3)(3A-1)V2^^^^ 1^ v/2(3A^ ) 

^ + ^ 2/n(-(3A-l)-(A + l)/(t))- 



(3.27) 



with 



_4(A-3)(3A- 1 

- 4(;,_3)(3;,_ 1)3/2 ^^ (3 + fim + m) - A(3 + mf) 

where 



_ M4t(l-3A)) 

4(3A-1) ^^-^^^ 

and C2 being the constant of integration. 

At this point we observe, that the above solution depends only on A. None of the parameters a 
or C, appear in equation (13.251) . As it seems the declination from general relativity in the kinetic 
term, produces more dramatic changes in the results. In the case where we restored the original 
kinetic term with A = 1 we got a solution which gives the classic one for type III when the other 
parameters assume the right values. In the latter case we turned off the extra potential terms and 
kept A as it is, resulting to a totally different solution. 



4 Quantum case for axisymmetric type III 



The Bianchi type III is a class B model. For these models it is known that the Euler-Lagrange 
equations from the Einstein - Hilbert reduced action with the scale factors taken as the degrees of 
freedom, are in general not equivalent to the reduced Einstein's equations. Specifically for type III 
this problem is solved in the diagonal case with ^ai3 = diag{a^, b"^, c^) by taking the axisymmetric 
condition a = c (dictated by the linear constraint). So, a first task is to check, if the same holds 
for the reduced Lagrangian emanating from (12.11) and the corresponding equations (13. 2p . (13. 3p . 
f)3.4p and (13. 5p . If we insert the scale factor matrix 

/a(t)2 
7a/3 = b{tf 




A 
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into action fl2.ip . we are led to the reduced Lagrangian 

bci^ CL^b'^ CLQjh Nh 

L{a, b, N, d, h) = 2a(l - 2A)— + a(l - A)—- - 4aA— - 2CNb + 2(C + 2r])— + aNa%{4.1) 

i\ bl\ i\ a-^ 

As it can be straightforwardly verified, the above Lagrangian gives rise to Euler - Lagrange equa- 
tions equivalent to the equations of motion (13. 2p . fl3.3p =( l3^ and fl3.5p (translated of course by 
Uiit) = 0, U2{t) = 21n(^), %(t) = ln(a(t))). We now proceed to the Hamiltonian formulation 
and the subsequent canonical quantization of this system. The conjugate momenta are given as 

Hat = (4.2) 
_ 4a((2A-l)&d-Aa&) 

2a{{X - l)a% - 2Xabd) 

= m ^^-^^ 

and thus the IIjv ~ is the primary constraint of the system. The Dirac - Bergmann algorithm for 
constrained systems ([20], [21]) leads to the canonical Hamiltonian, He = alia + bUb — L = NT-Lc, 
with 

n - + j2A-2)^ X 2{C + 2r,)-2ia^ + aa^ 

^^"8a(3A-l)b^^" + 4a(3A-l)a2^^^ 2a(3A - l)a ^'^^^^ ^ ^^'^^ 

The requirement of preservation in time of the primary constraint indicates that "He ~ is the 
secondary constraint. The algorithm thus terminates for A 7^ | and we are left with these two first 
class constraints. The case A = | will be treated separately. One can read off the supermetric 
from the kinetic part of He- 

( A \ 

_ 4a(3A-l)b 2a(3A-l)a 
^ ~ A (2A-1)6 



\ 2a(3A - l)a 2a(3A - \)a^ ) 
To canonically quantize the system in the Schrodinger picture we adopt the usual choice of basic 
operators [fi = 1): Iliv Hat = Hq — )• Ha = lib ^ t^b = with associated 

commutation relations [A^, 11'^] = iSnn' \cL^Yl'g\ = iSaa', [^5 n'fe] = i^bb'- 

We then follow Dirac's proposal and define the wave function to be annihilated by the quantum 
analogue of the two first class constraints. The primary n7v^(A^, a,b) = informs as that does 
not depend on A^. For the secondary we choose the factor ordering so that the kinetic part of Tic 
becomes the Laplace - Beltrami operator 

nk^n = ^r-i/^n^ri/^r'^-n^ (4.6) 

The Wheeler - DeWitt equation T-Lc^i^a, 6) = becomes 

A-1 (2A-1)6 d^^ X d^^ , 

H V (4.7) 



8a(3A - 1)6 da^ 4a(3A - \)a^ dh^ 2a(3A - l)a dadb 
A-1 ^ 2A-1 ^2{C + 2r])-2ia^ + aa^^_^ 



8a(3A - l)ab da 4a(3A - l)o? db a? 
The space of solutions to this partial differential equation is difficult enough to be found in full 
generality (i.e. for arbitrary values of the parameters). Thus, we give below explicit solutions for 
A = 1, cr = and A = |. 
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4.1 Case A = 1 

Under the assumption A = 1 (from now on we take a = 1) equation (14. 7p reduces to 

b 1 d^^ ^ 1 <9* ^ 2(C + 277) -2^a2 + aa^^_ ^ ^^^^ 



Under the change of variables (a, b) i— t- (x, y) where x = ab and ?/ = we can write the previous 
equation without a mixed derivative term. A further simphfication by dropping the cosmological 
term, i.e. setting a = 0, results in the Wheeler - DeWitt equation being solvable through separation 
of variables: \E'(x,?/) = X{x)Y{y) whith the two one variable functions satisfy the equations 

x'^X"(x) xX'(x) ^ o o , , 

y'Y%y) , yY\y) , 2(C + 2r^) 

8Y{y) + 8Yiy) + y- ' ^^'^^^ 

with c being an arbitrary constant and the primes indicating differentiation with respect to the 
arguments. The solutions to these equations are 

X{x) = HiM^^/^x) + K^Y^A^/^x) (4.12) 
Y{y) = - 4c)/,4c(4 ^~^ ~ ^ -^) + K4(-1)^T(1 + 4c)/4e(4^^^^^^) (4.13) 

y y 

where Kj, i = 1,2,3,4 are the constants of integration Jy(x), Y^{x) and Iu{x) are the Bessel J, 
Y and I functions respectively, and lastly T{x) stands for the Gamma function. Translating this 
result to the old coordinates the solution of the Wheeler - DeWitt without a and for A = 1 is 

^(a,6) = (/ciJ4c(4v^a6) + /C2l4c(4v^a6)) («:3(-l)"'T(l - 4c)/_4e(4V-C - 2t]) + 

+A€4(-1)'T(1 + 4c)/4c(46v/-C - 2r/)) (4.14) 



If we want to compare with the results from Einstein's gravity, we must start from the reduced 
Lagrangian 

2ba? Aaab /. , 

L, = -^-^-2Nb iA.15) 

We follow exactly the same procedure, which results in the following Wheeler - DeWitt equation: 



8a2 Aadadb %a? db 

Again, the choice of variables x = ab and y = leads to the extinction of the mixed derivative 
and allows fl4.16p to be solved with a separation of variables. Upon setting '^{x,y) = X{x)Y{y) 
we get: 

x^XUx) ^ xJ^x) ^ ^ 2c2 (4.17) 



8Xe(x) 8Xe(x 



+ "f,: = 2c' (4.18) 



8YE{y) 8YE{y) 
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These are the same equations one would acquire by setting the values (^ = l,C = ^ = Oor 
^ = 1, C = -2ri in (H:TUD - (HTTD . The solution of (HTTBD is: 

^^(a, 6) = (/€iJ4c(4a6) + K2>4c(4a&)) ^/ta cosh(4cln(^)) +i K4 sinh(4cln(^))^ (4.19) 

The first parenthesis is the same with the one from the first solution under the condition = 1. 
The terms in the second parenthesis are now written in the form of ordinary functions, as opposed 
to the previous case. If we now set ,^ = 1, C = ?7 = in (14.141) . we get 



{^2 + K3) (/ti Jo(4a6) + K2Yo{Aab)) if c = 
60 if |4c| G M+ 



^{a,b) = 

while at the same time for c — 

(c^o)(a, b) = Ks (KiJoiAab) + K2Yo{Aab)) (4.20) 

We have to note here that fl4.20p is not produced by the solution of the set of equations (14.171) 
and (14.181) with c = 0. We just refer to the limit of their solutions with c 7^ and as c — ?■ 0. So 
under the limiting case, where the extra coupling constants tend to take values that correspond 
to Einstein's gravity, we see that "^{a, b) diverges, unless we make the specific choice in (I4.10p and 
(14. lip to set c = 0. Then the solution becomes essentially the same with limc_>o ^£;(a, &), which 
of course is not the general solution of (I4.17P and (14.181) for c = 0. The latter is e ic=o){(i,b) = 
{ks + K4 In i) (niJoiAab) + K2^o(4a6)). 

One could remove the divergence of \I'(a, b) as C and rj tend to zero and c 7^ 0, by extinguishing 
one of the constants of integration K3 or K4 (the choice depends on whether the admissible values 
of c are positive or negative), but that would eventually lead to a zero wave function. 

4.2 Case A = | 

Equation (14. 5 p entails the need to treat the case A = | separately. If we set this value for A in 
(14. ip and calculate the conjugate momenta Ua and Ilf, we see that they are not independent, a 
fact which results in the existence of another constraint. It is easy to check that in our case 



or re-writing it in operator form 



aUa + bUb = (4.21) 



An invariant function for this operator is ^, so if we set this fraction as a new variable, the operator 
will be cast into canonical form, and the transformed Lagrangian will depend on only one velocity. 
We choose the set of new variables to be (u, v) with u = ^ and v = a. Under this transformation 
and for the specific value A = | Lagrangian (14. ip becomes 

2ci'^'ii^ / u \ 

Liu, a, A^, ii) = — — + N (2(( + 2r])- - 2^ua + aua^ (4.23) 
6Nu \ a / 
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The conjugate momenta are 

Hiv = (4.24) 

Ha = (4.25) 
^ 4a^-u 

n„ = — — 4.26 

3Nu ^ ' 

with Hat and Ha being the primary constraints. Now, we can write the canonical Hamihonian 
He = iiHu — L, which turns out to be 

= NTic = Nu i^Iil - 2i^±M + 2ia - aa^\ (4.27) 
\6a-^ a J 

while the total Hamiltonian defined in the full phase space is 

Ht = H, + ^iHiv + li2Ha (4.28) 

The condition that the constraints must be preserved in time, {Y\.n,Ht} ~ and {YlaiHT} ~ 
leads to the secondary constraints 

= 16{C + 2ri)a^ - W^a^ + 8aa^ - 3nl (4.29) 
X2 = 2(C + 2?7) -4a2^ + 3aV (4.30) 

As for their preservation in time, the Poisson bracket of Xi with the total Hamiltonian is weakly 
zero, while {x2,Hx} ~ imposes the condition /i2 ~ in 04.281) . So we have a total of four 
constraints with Hat, xi being first class and Ha, X2 being second class. We use the definition 
of the Dirac bracket { , }d = { , } ~ { } with being the elements of the set of 

the second class constraints and A^j, = the non singular matrix formed by their Poisson 

brackets. 

The use of the Dirac brackets allows us to set the second class constraints strongly equal to zero. 
The basic canonical commutation relations become 

{N,I{n]d = 1 

{w,n„}z) = 1 

{a,njz) = 

with the last of them simply emphasizing the redundancy of the corresponding pair of variables. 
Thus on the further restricted phase space X2 — ^ leads to a constant scale factor 

« 2e±V4e-6(C + 2,). 

3(7 

It is interesting to note that, as expected from the equivalence between the Hamiltonian and the 
Lagrangian formulation for constrained dynamics, (14.311) is identical to the value of the square 
of the constant scale factor in the classical solution, as obtained in (13.221) . Again we choose the 
operators 

^ ^ d 

H„ H„ = -I— 
ou 
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The quantum analogue of the primary constraint nAr\l/(A^, m) = dictates a solution of the form 
\E' = \E'(m). Our Wheeler - DeWitt equation results from xi\E'(w) = under the condition fl4.3ip . 

32 (a^^ - 9(C + 2r])Ca ± 2^^^/AC'^ - 6{C + 2r])a T 3(C + 2ri)a^/4:C'^ - 6{C + 2r])a^ - 81a^^" 

its solution is 

^{u) = cie^^" + cae"^" (4.33) 

with 

32 (4^ - 9(C + 2r,)(a ± 2e ^W^W+M^ T H( + 2r,}a^ie " 6(C + 2)))ff) 

81^^ f^-^^) 

Upon setting ( = —2r] and ^ = 1, for the plus value of f l4.3ip . we are led to k = On the other 
hand, the minus value gives k = 0, which induces a constant wave function for every value of u. 



u 



5 Discussion 

We have treated a diagonal Bianchi type III cosmological model within Horava - Lifshitz theory 
classically and quantum mechanically. 

At the classical level, the main tools for investigating the solution space were: (a) the automor- 
phisms of the corresponding Lie algebra whose generators constitute Lie point symmetries of the 
equations of motion and (b) the use of the quadratic constraint as an algebraic equation deter- 
mining the lapse N{t) and the subsequent replacement of its value in the spatial equations. The 
virtue of this is twofold: On the one hand it becomes possible to reduce the order of the system of 
differential equations and ultimately arrive at f l3.6p and the main solution fl3.14p . and on the other 
hand various signature solutions are attained. When the lapse is not determined by the quadratic 
constraint, we are led to the marginal metrics f l3.10p and fl3.19p . These metrics are interesting 
enough, since they are curvature homogeneous spacetimes. Additionally, one might consider fl3.19p 
for (Lt2 = as equivalent to fl3.1Up . since the curvature invariant relations are then identical. But 
this is not true, because the transformation needed for this identification involves mixing of t,z 
as it can be seen by the corresponding killing vectors. Thus the metrics are not equivalent in the 
context of Horava - Lifshitz theory. For the same reason line element fl3.19p can not be considered 
as static. It is also interesting that for A = | the resulting line element fl3.24p is also a curvature 
homogeneous space despite the fact that in this case the lapse has been determined. For these 
three spacetimes a noteworthy segregation of the original coupling constants occurs: if we change 
their values in a way that does not affect the combinations entering the essential constants the 
classical solutions are not affected, leading one to consider all these models as equivalent. 
Another important observation concerns the coupling constant A. Its value significantly changes 
the essence of the classical theory. For A = 1, the general solution f l3.14p can be considered as 
the original solution of Einstein's equation plus higher term "corrections" . For A = | the solution 
(I3.24P develops for cr = curvature singularities at the Einsteinian limit, or becomes a solution of 
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Einstein's equation plus a cosmological constant term. The only case for which A becomes utterly 
unimportant is the static pseudometric f l3.10p . a thing expected since A appears inside the kinetic 
term. Lastly, it is fortunate that, even though we consider a class B model for which we can write 
a valid Lagrangian only under the axisymmetric condition, we do not lose this good property for 
the action of Horavas' theory. Thus the Euler - Lagrange equations of the reduced action, are 
identical to the reduced equations of motion. 

At the quantum level, we were able to solve the Wheeler - DeWitt equation for the values A = 1 
and A = |. Even in the case of A = 1 the Einsteinian limit of the wave function can not reproduce 
the corresponding Einsteinian wave function, because it diverges. The only non divergent case is 
when the separation constant is set to zero. 

Of course, the renouncement of four dimensional covariance is not a problem free choice: it has 
been observed ([5], [10], [22]) that there are not only difficulties attaining general relativity in the 
Einsteinian limit (A — )■ 1) but also stability problems regarding the UV region. Recently, a new 
approach is proposed in [2^, in an attempt for the theory to overcome its problems. Since we 
could not obtain the general solution of the Abel equation (13. 6p . we can not comment on this 
issue. 
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